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We consider a linear O(N) model and introduce an auxiliary field to eliminate the scalar self- 
interaction. Using a suitable limiting process this model can be continuously transformed into the 
nonlinear O(N) model. We study this model at nonzero temperature applying the CJT formalism. 
Due to the presence of the auxiliary field the effective potential and, consequently, the equations for 
the order parameter and the masses of the particles turn out to be different from those previously 
obtained in similar studies of the linear and nonlinear O(N) models. The order of the chiral phase 
transition depends sensitively on the choice of the renormalization scheme. In the linear version of 
04 ■ the model and for explicitly broken chiral symmetry, it turns from crossover to first order as the 

mass of the a particle increases. In the nonlinear case, the order of the phase transition turns out to 
be of first order. In the region where the parameter space of the model allows for physical solutions 
of the gap equations, Goldstone's theorem is always fulfilled. 
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I. INTRODUCTION 



O i 1 Scalar models in d + 1 space-time dimensions with orthogonal symmetry are widely used in many areas of physics. 
JL | Some applications of these O(N) models are quantum dots, high-temperature superconductivity, low-dimensional 
systems, polymers, organic metals, biological molecular arrays, and chains. In this paper, we focus on a physical 
| system consisting of interacting pions and sigma mesons at nonzero temperature T. For three spatial dimensions, 
l—— d = 3, an analytical solution to this model does not exist. Thus, one has to use many-body approximation schemes 
in order to compute quantities of interest, such as the effective potential, the order parameter, and the masses of the 
particles as a function of T . As an approximation scheme never gives the exact solution, it is of interest to compare 
00 1 different schemes and assess their physical relevance. 

t-H ! For N — 4 the O(N) symmetry group for the internal degrees of freedom is locally isomorphic to the chiral 
OS ■ SU{2)r x SU(2)l symmetry group of quantum chromo dynamics (QCD) with Nf = 2 massless quark flavors. The 
[ phenomena of low-energy QCD are largely governed by chiral symmetry. 

In the case of zero quark masses the QCD Lagrangian is invariant under U(Nf)n x U{N})l transformations, Nf 
being the number of quark flavors. However, the true symmetry of QCD is only U(Nf)v x SU(N/)a, because of the 
axial anomaly which explicitly breaks U(1)a due nontrivial topological effects For Nf nonzero but degenerate 
quark masses, the SU(Nf)A symmetry is explicitly broken, such that QCD has only a U(Nf)v flavor symmetry. In 
J> , reality, different quark flavors have different masses, reducing the symmetry of QCD to U(l)v, which corresponds 
? ■ to baryon number conservation. In the vacuum, the axial SU (Nf)A symmetry is also spontaneously broken by a 
non- vanishing expectation value of the quark condensate (qq) ^ 2] . According to Goldstone's theorem, this leads 
■ to Nf — 1 Goldstone bosons. 

The chiral symmetry is restored at a temperature T which for dimensional reasons is expected to be of the order 
of Aqcd ~ 200 MeV. This scenario is indeed confirmed by lattice simulations, in which a crossover phase transition 
at T c ~ 150 MeV has been observed @. 

For vanishing quark masses, the high- and the low-temperature phases of QCD have different symmetries, and 
therefore must be separated by a phase transition. The order of this chiral phase transition is determined by the 
global symmetry of the QCD Lagrangian; for U{Nf)y x U(Nf)A> the transition is of first order if Nf > 2; for 
U{Nf)v x SU(Nf)A the transition can be of second order if Nf < 2 4]. If the quark masses are nonzero, the 
second-order phase transition becomes crossover. 

The calculation of hadronic properties at nonzero temperature faces serious technical difficulties. For a nonconvex 
effective potential standard perturbation theory cannot be applied. Furthermore, at nonzero temperature an additional 
scale is introduced which invalidates the usual power counting in terms of the coupling constant |5|. A consistent 
calculation to a given order in the coupling constant then may require a resummation of whole classes of diagrams 
i- 

A convenient technique to perform such a resummation and thus arrive at a particular many-body approximation 
scheme is the so-called two-particle irreducible (2PI) or Cornwall- Jackiw-Tomboulis (CJT) formalism [3, Q, which 
is a relativistic generalization of the ^-functional formalism [t| [l(| ■ The CJT formalism extends the concept of 



2 



the generating functional T [(f)] for one-particle irreducible (1PI) Green's functions to that for 2PI Green's functions 
r [</), G], where <j) and G are the one- and two-point functions. The central quantity in this formalism is the sum of all 
2PI vacuum diagrams, T 2 [4>, G]. Any many-body approximation scheme can be derived as a particular truncation of 

r 2 [0,G]. 

An advantage of the CJT formalism is that it avoids double counting and fulfills detailed balance relations and 
thus is thermodynamically consistent. Another advantage is that the Noether currents are conserved for an arbitrary 
truncation of T2, as long as the one- and two-point functions transform as rank-1 and rank-2 tensors. A disadvantage 
is that Ward-Takahashi identities for higher-order vertex functions are no longer fulfilled [ll|. As a consequence, 
Goldstone's theorem is violated [H, EH- A strategy to restore Goldstone's theorem is to perform a so-called "exter- 
nal" resummation of random-phase approximation diagrams with internal lines given by the full propagators of the 
approximation used in the CJT formalism [TTj . 

In the literature different many-body approximations have been applied to examine the thermodynamical behavior 
of the O(N) model in its linear and nonlinear version. In Ref. [l4j optimized perturbation theory was used to compute 
the effective potential, spectral functions, and dilepton emission rates. The CJT formalism has been applied to study 
the thermodynamics of the O(N) model in the so-called "double-bubble" approximation [13, 0, , in Ref. (24{ 

sunset-type diagrams have been included. The I/TV expansion has also been used several times to study various 
properties of the O(N) model at zero [III [26| and nonzero (27l - [3Cl | temperature. 

Here we calculate the effective potential, the masses, the condensate, and the pressure of the O(N) model at nonzero 
T both in the linear and nonlinear cases by making use of an auxiliary field. This field allows us to obtain the nonlinear 
version of the model by a well-defined limiting process from the linear version. We find that to lowest loop order, 
T2 = 0. Nevertheless, when using the gap equation for the auxiliary field, the gap equations for pions and sigma 
mesons still contain self-consistently computed loops. Although qualitatively similar to the standard double-bubble 
approximation in the treatment without auxiliary field, the gap equations are quantitatively different and thus lead 
to different results for the order parameter and the masses of the particles as a function of T. 

The order of the chiral phase transition depends sensitively on the choice of renormalization scheme. In the linear 
version of the model and for explicitly broken chiral symmetry, it turns from crossover to first order as the mass of 
the a particle increases. In the counter-term renormalization scheme, this transition happens for smaller values of 
the a meson than in the so-called trivial renormalization scheme. In the nonlinear case the phase transition is of 
first order. Besides, in the region where the parameter space of the model allows for physical solutions of the gap 
equations, Goldstone's theorem is always respected. 

The manuscript is organized as follows: in Sec. II the linear and nonlinear versions of the model are presented and 
how they can be related with the help of an auxiliary field. In Sec. Ill the effective potential and the gap equations 
are derived. In Sec. IV the results are presented for the linear and nonlinear versions of the model in the case of 
non-vanishing and vanishing explicit symmetry breaking. Section V concludes this paper with a summary of our 
results and an outlook for further studies. 

We use the imaginary-time formalism to compute quantities at nonzero temperature. Our notation is 

f(k)=T / —3/(2™^) . (1) 

We use units h = c = fcg = 1- The metric tensor is g^ v — diag(I, —I, —I, —I). 



II. THE O(N) MODEL 
The generating functional of the cr-model with O(N) symmetry at nonzero temperature T is given by 

rl/T 



Z L (e,h) — N JvaV^ exp ^ 



dr / d 6 xC^ a , (2) 



v 



with the Lagrangian 



C a - a = \{d^f -U($,a) , U($,a) = ^a($ 2 -v 2 ) + ^a 2 -h* , (3) 

Z Z o 

where <& 2 = $* = (<7,7Tl, . . . , ttn-i) and a is an auxiliary field serving as a Lagrange multiplier. One can obtain 
the generating functional of the O(N) model in its familiar form by integrating out the field a: 



Z L (e,h) = J V<S> exp uf 



l/T 

dr I d 3 xC a I , (4) 
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with the Lagrangian 



<S> 2 -vl) 2 



ha 



(5) 



As one can see, the potential of the model exhibits the typical tilted Mexican hat shape, with the parameter 1/e being 
the coupling constant, h the parameter for explicit symmetry breaking, and vq the vacuum expectation value (v.e.v.) 
of $. The tt{ holds can be identified as the pseudo-Goldstone fluctuations. 

Another way to see the equivalence to the standard form of the O(N) model is to use the equation of motion for 
the auxiliary field a, 



6a 



0, 



6£<0-OL 

Sd^a 







Ne 



{&-v 2 ) 



(6) 



When plugging the latter into £ CT _ Q one recovers, as expected, the familiar Lagrangian C a . 

The advantage of the representation @ of the generating functional of the linear a- model is that, by taking the 
limit e — > + , one naturally obtains the nonlinear version of the a- model with the fields constrained by the condition 

$2 



Vq. In fact 





' f 1/T f 


J Va X>$ exp 


/ dr d 3 x£ a . a 




Jo Jv 



Jv<f>5{<P 2 -v 2 )exp ' dT J d 3 x Q (<9 M $) 2 + ha 



(7) 



because (5(<E> 2 — Vq) can be identified with 







j T>a exp < 









(8) 



which is the mathematically well-defined (i.e., convergent) form of the usual representation of the functional 
5— function. Equation ([5) ensures that the Mexican hat potential becomes infinitely steep and, consequently, the 
mass of the radial degree of freedom infinite. 

Note that in some previous studies of the O(N) nonlinear a— model [27l. [28j. the e-dependence in Eq. ((5J) was not 
appropriately handled: there, the limit e — > + was exchanged with the functional T>a integration, effectively setting 
e = in the exponent. This, however, is incorrect, since the additional term ~ e a 2 is essential to establish the link 
between the linear model and the nonlinear one. Without this term, an integration over the auxiliary field does not 
give the correct potential of the linear model. Thus, for a proper construction of the nonlinear limit of the O(N) 
model the e-dependence must be included. 



III. THE EFFECTIVE POTENTIAL AND GAP EQUATIONS 



In this work we study the thermodynamical behavior of the O(N) c-model, and in particular the temperature 
dependence of the masses of the modes and of the condensate. To this end one has to apply methods that go beyond 
the standard loop expansion which is not applicable when the effective potential is not convex [3l|, as is the case 
here because of spontaneous chiral symmetry breaking. A method that allows to compute quantities like the effective 
potential, the masses, and the order parameter at nonzero temperature is provided by the Cornwall- Jackiw-Tomboulis 
(CJT) formalism [tJ. Within this formalism, the effective potential V is given by 



V cS = U{<j ) ) + ^J[\nG- 1 {k) + D- 1 {k-A)G(k)-l\ + V 2 {^G) , 



(9) 



where the one-point function <fi and the connected two-point function G(k) are determined by the extrema of the 
effective potential, 



sv eB 







5V C 



cir 



6G(k) 



. 



(10) 
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In Eq. (j9)), U(4>) is the tree-level potential, D~ 1 (k;4>) is the tree-level propagator in momentum space, and Vi{<(>,G) 
contains all two-particle irreducible diagrams. 

In our case, the fields occurring in the Lagrangian are a, 7r = (iri, . . . , 7Tjv-i)* , as well as the auxiliary field a. The 
tree-level potential is given by 

1 Ne 

U = 7r(ao + a) (a 2 + tt 2 + 2(70 + 2 - v 2 ) + —(a + a) 2 - h{4> + a) , (11) 

2 8 

where we have already anticipated that the fields a and a attain non- vanishing vacuum expectation values, a — > cf> + <7 
and a — > ao + a, respectively. The shift produces a bilinear mixing term, iacr<f>, which renders the mass matrix 
non-diagonal in the fields a and a. There are various ways to treat this mixing term, for instance keeping this term 
and allowing for a non-diagonal propagator which mutually transforms the fields a and a into each other. There is, 
however, a simple method to eliminate this unphysical mixing: we perform a further shift, 

i<t> 

a — >a-A-±-a. (12 

Ne v ; 

Using the fact that linear terms in the fields vanish on account of the definition of the vacuum expectation values, 
dU /d(f> = 0, dll/dao = 0, the resulting expression of the Lagrangian reads 

r 1 m ^ 1,~ ,2 1/, 40 2 \ 2 I 2 1 Ne 2 

£«- a = 2 ( m^) + 2 (^ 7r ) - 2 \ la ° iVe J °" ~ 2 ^ 1 

_ + tt 2 ) - a(a 2 + tt 2 ) - ^a o (0 2 - ug) - ^a 2 + H ■ (13) 

One can immediately deduce the inverse tree-level propagators and the tree-level masses: 

Dr 1 (k;<j),ao) = -k 2 + m 2 ; i = cr , tt , (14) 
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40 

ia ° + ~Ne ' TO ^ = ia ° ' 



(15) 



_/V"e 

D- 1 (k;^a ) = m 2 a = — . (16) 

The shift (|12p has two features: the Jacobian associated with the transformation is unity, thus the functional inte- 
gration in Eq. ([7]) remains unaffected. Furthermore, this shift generates a term in the c-mass, which diverges in the 
limit e — > + , see Eq. (|15|) . This is expected, since the u-particle becomes infinitely heavy in the nonlinear version of 
the 0(N) model. 

We now turn to the effective potential, 

V cS (<l>,a ,G a ,G„,G a ) = ^a o (0 2 - v 2 ) + ^a 2 - lup (17) 



\ £ /[lnGr 1 (fc)+i?r 1 (*;^«o)G < (fc)-l]+% 



2 

We now make an important observation. The Lagrangian (|13|) only contains three-point interaction terms between 
the fields a, re, a. Therefore, the lowest-order contribution to V% is a two-loop "sunset"-type diagram constructed 
from two three-point vertices. There is no two-loop diagram of double-bubble-type which would require the existence 
of a four-point vertex. If we restrict ourselves to a many-body approximation where the self-energy of the particles 
is independent of momentum, we have to neglect sunset-type diagrams (as well diagrams of higher loop order), and 
the contribution of V% to the CJT effective potential vanishes trivially. As a consequence, the effective potential only 
contains tree-level and one-loop terms. 

Using the stationary conditions for the effective potential 

SV c g fiVeff dVcff 

one derives the following equations for the condensates, 

h = ia <j)+^ f G a {k) , (19) 



k 



Ne 



-v 2 + G a (k) + (N-1) / G v {k) 



(20) 
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and full propagators, 



G-\k) = -k 2 + M? = Dr\k; 0, a ) , 



Mf = m 2 ; 



i = a , 7r , a . 



(21) 



Since V% vanishes in our approximation, the contribution from the 1PI self-energy, If; = — IdVij dGi, trivially vanishes 
and the full propagators are equal to the tree-level propagators. The thermodynamic pressure is, up to a sign, identical 
to the value of the effective potential at its minimum, 

p = -V e s((f>,ao,G a ,G„,G a ) , (22) 

where (j), cto, Gi are the solutions of the gap equations (fT9]). (|20|). (|2"T]) . 

The effective potential in Eq. ()1T[) depends on the physical field 4> and the auxiliary field ao. Furthermore, it has the 
shape of a saddle. Therefore the solutions of the gap equations (IT9)) . (|2T))) . (|2"Tj) correspond to a saddle point and not 
to the global minimum of Eq. (117[) . Since the auxiliary field is not an independent dynamical degree of freedom but a 
Lagrange multiplier, the condensate equation for ao serves only as an additional constraint on the thermodynamics 
of the model. Substituting iao in Eqs. (fT5)) and (flU| by the right-hand side of Eq. (f2T)|) gives the following equations 
for the physical condensate <j) and the masses 



M 2 (e,h) + 



M 2 (e,h) = M 2 (e,h) + 



Ne 

2 



G a {k) 



Ne 



v 2 + / G a (k) + (N-1) / G«(k) 



(23) 



This set of equations will be referred to as "Hartree approximation" in the following. Neglecting terms which are 
subleading in 1 /N — an approximation commonly referred to as the large- N limit — the gap equations reduce to 

h = (^M 2 (e,h) + 0(N- 1 ) , 



M 2 (e,h) = M 2 (e 1 h) + ^- , 



Ml{e,h) 



Ne 



N 



Note that the condensate and the v.e.v. are 



N, i.e., 



- N. 



OiN- 1 ) 



(24) 



Furthermore, after substituting the auxiliary field in Eq. (fT7f using Eq. (|20|) . we obtain the usual "Mexican hat" 
shape for the effective potential: 



V cS {<t>,G ai G^) 



1 



2Ne 
1 

2Ne 



G a {k) 



hcj) 



2 

(n-K 

2Ne 



'I — (T,7T 

2 r 



l\nG^(k) + Dr\k;<j))G t (k) - l] 
N - 1 



Ne 



GAk) / G a (q) . 



(25) 



Finally, we present for comparison the results for the effective potential that one obtains by a direct application 
of the CJT formalism to the original Lagrangian C a of the O(N) linear cr-model, Eq. ([5]), where no auxiliary field is 
present: 



K£ JT (0,G ff ,Or) 



1 



2Ne 



,2\2 



^[lnGr 1 (fc) + J Dr 1 (A : ;0)G l (A : )-l] + K i CJT (G ff ,G,) , 



(26) 



where 



V^' J1 (G ff ,G.) 



2Ne 



G a {k) 



+ (iV + l) 



N- 1 

2Ne 



G„{k) 



TV- 1 
Ne 



G„(k) / G a (q) . 



(27) 



The expression for Vp is rather similar to the last line of Eq. ([25]) . except for the combinatorial prefactors. However, 
as we shall comment in the next section, the numerical results using the effective potential (1261) differ sizably from 
the ones obtained using the effective potential (|2"5"|) . 
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In Ref. [32], Ivanov, Riek, and Knoll introduced a heuristic modification of the CJT result (f25j) . which is constructed 
such that Goldstone's theorem is preserved. In our notation, 



V 2 (Go-jGtt) — V 2 (GujG-k) 
5-2N " 



N — 1 



2Ne 



Ne 

2 



[G„(k) - G a {k)} 



G a {k) 



{N-lf 
2Ne 



G„(k) 



3(N-1) 

Ne 



G„(k) / G a (q) 



(28) 



Note that the coefficient of the term containing two pion loops is the same as in Eq. (|25p. although all other terms have 
different coefficients. However, in the large- N limit, all approaches give the same result for the effective potential. 



Carrying out the Matsubara summation of the tadpole integrals gives 



Gi(k) 



d 3 k 



1 



2\/h? + M? 



exp 



k 2 + Mf/T - 1 



(29) 



Note that the vacuum contribution of Eq. (|29p can be written as 

f d 3 k 1 f d 4 fc 



{^f 2y/P + M 2 J (2nfk 2 +M 2 ■ 
This term, however, exhibits quadratic and logarithmic divergences and thus has to be regularized, 



j£k_ 
(2tt) 4 



M 2 



k 2 + ^ 2 



f-M 2 
(fc 2 + /x 2 ) 2 



(4tt) 5 



M 2 {n (M 2 
M 2 V V 2 



M 2 



(30) 



(31) 



where fi is the renormalization scale. We choose to put fx = m a for the following reason. In the vacuum and in the 
chiral limit, h — > + , the pion mass vanishes, —¥ + . Then, the first equation (|23[) requires that the vacuum 
contribution to the tadpole integral J, G a (k) vanishes. From the result (|3"Tj) we observe that this can only be realized 
with the choice /i = m a - Note that is always positive semi-definite, irrespective of the value of M 2 / ji 2 . 



IV. RESULTS 



In this section, we show numerical solutions of Eqs. (|23|) and (|24|) for N = 4, corresponding to a system of three 
pions and their chiral partner, the a particle. We discuss the results for the linear and the nonlinear ct- model, with 
and without explicitly broken chiral symmetry. Furthermore, we investigate different renormalization schemes, i.e., 
the above mentioned counter-term renormalisation (CTR) method, Eq. (|3Tj) . and a so-called trivial renormalisation 
(TR) where the vacuum contribution Q M of the tadpole integral is set to zero. 



A. Linear model with explicitly broken symmetry 

In Fig. [T] we show the masses of the pion and the sigma, as well as the condensate as a function of temperature 
for different values of the vacuum sigma mass m a . One observes that the condensate decreases as a function of 
temperature, which is a consequence of the restoration of chiral symmetry. Depending on the value of m a , chiral 
symmetry restoration may proceed via a phase transition. In the CTR scheme, the phase transition is of second order 
for m a ~ 500 MeV, and of first order for larger values of m a . For smaller values, however, the transition is only 
crossover. In the chirally restored phase, the condensate is always nonzero because of the small explicit breaking of 
chiral symmetry due to non- vanishing quark masses (or, equivalently, m, = 139.5 MeV). Since the results for the 
TR scheme are qualitatively similar, we do not show them explicitly, but we remark that the second-order transition 
occurs for larger values of the vacuum sigma mass, m a ~ 700 MeV. Note that a crossover transition is also found in 
lattice QCD calculations Q. However, it should be stressed that in QCD the identification of the chiral partner of 
the pion is a long-debated issue, see Refs. [33[ and refs. therein. 

Figure [5] shows the effect of different renormalization and approximation schemes on the behavior of the masses and 
the condensate as functions of temperature. We keep the vacuum mass of the sigma fixed to m a = 550 MeV. In the 
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FIG. 1: The pion mass, the sigma mass, and the condensate as a function of T in the 0(4) linear model in case of explicitly 
broken symmetry using the CTR scheme for different values of m a . 




FIG. 2: The pion mass, the sigma mass, and the condensate as a function of T in the 0(4) linear model in case of explicitly 
broken symmetry for m CT = 550 MeV and different renormalization schemes. 



CTR scheme, the system exhibits a first-order phase transition. When using the TR scheme, however, one observes 
a crossover transition. In the large- N limit with CTR, the chiral transition is always crossover, independent of the 
mass of the a. In Fig. [2] the crossover transition is observed to be smoother for the large- N approximation with CTR 
than for the other cases. In fact, with this renormalization scheme, the smoothness is proportional to m a . We shall 
see in the next section that the transition disappears as we approach the nonlinear limit m a — > oo. This, however, 
does not happen for the TR scheme. 



B. Nonlinear model with explicitly broken symmetry 

The results for the nonlinear model are obtained by taking the limit e — s- + . Because of the relation (|15l) . 1/e = 
(m^. — mj) /4> 2 , the nonlinear limit is equivalent to sending m a to infinity. 

In this case, in the Hartree approximation and when the TR scheme is used, the phase transition is of first order, 
with a rather large discontinuity in the condensate at the critical temperature of T c ~ 178.6 MeV, see Fig. [3] The 
condensate is very small above T c , but still nonzero, because of explicit symmetry breaking. The first-order nature of 
the transition is in line with the expectation from the linear case, where the transition becomes first order when the 
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FIG. 3: The pion mass and the condensate as a function of T in the 0(4) nonlinear model in case of explicitly broken symmetry 
using the Hartree-approximaion in the TR-scheme for m a oo (in practice m a = 250 MeV is used). The solid line corresponds 
to the physical case. The dashed line corresponds to the unstable or metastable solution of the gap equations. 




50 100 150 200 250 300 350 400 50 100 150 200 250 300 350 400 

T [MeV] T [MeV] 

FIG. 4: The pion mass and the condensate as a function of T in the 0(4) nonlinear model in case of explicitly broken symmetry 
using the LN-approximaion in the TR-scheme (full) and CTR-scheme (dashed) for m CT — s> oo. 



a mass is sufficiently large. Below T c the a mass is infinitely heavy and there are only pionic excitations. Above T c 
the masses of a and 7r become degenerate. 

In the large- TV approximation and with the TR scheme, the phase transition is crossover with T c ~ 185 MeV, see 
Fig. |U In this case the a field remains infinitely heavy also above T c . This is the main difference to the Hartree case, 
where the a becomes degenerate with the tt above T c . It is at first sight surprising that this small difference can cause 
such a drastic change in the order of the phase transition. The explanation lies in a comparison of the gap equations 
(PI51) in the Hartree approximation with those in the large- N limit, Eqs. (f!?4")). Since the a is infinitely heavy below T c , 
there is no contribution from this mode to the gap equations. However, above T c , thermal fluctuations of the a can 
contribute in the Hartree approximation, while they remain absent in the large- N limit. This is sufficient to drive the 
transition to first order in the Hartree approximation. 

In the CTR scheme, the parameter space of the model does not allow for physical solutions in the nonlinear case 
m a — > oo. In the Hartree approximation, (f> — > and M ai — > oo for all values of T. In the large- N approximation, 
the transition disappears completely, the condensate and the masses retain their constant tree-level values for all 
T > 0: (j> = .U, M a = m ff , M„ = m„, see Fig. 0] 
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FIG. 5: The pion mass, the sigma mass and the condensate as a function of T in the 0(4) linear model in the chiral limit. 



C. Linear model in the chiral limit 



The chiral limit is obtained by taking — > + . In this case, the first equation (|23|) requires that the a— tadpole 
exactly cancels M%. This, however, is only possible, if the pion becomes tachyonic, M% < 0, since the thermal as well 
as the vacuum contribution (|3 1[) to the tadpole are always positive (semi-)definite. As a consequence, we can only 
show results in the large- N limit, since there this problem is absent, cf. Eqs. 

In Fig. [5] we show the behavior of the masses and the condensate as functions of temperature for various values of 
the vacuum sigma mass in the large- N limit in the CTR scheme (the results for TR scheme are qualitatively similar, 
therefore we do not show them explicitly). The results of Fig. [5] are in agreement with universality class arguments 
which predict a second-order phase transition. In the phase where chiral symmetry is spontaneously broken the pions 
are massless in accordance with Goldstone's theorem. Above the critical temperature the chiral partners become 
degenerate in mass. The condensate as a function of temperature is independent of the value of m a . This can be 
seen as follows. We subtract the last equation (|24|) at T = (where <fi — /„■) from the same equation at an arbitrary 
temperature T < T c , where T c is the phase transition temperature. Since in the phase of broken chiral symmetry we 
always have Mj = 0, the result is 

= tf(T)-tf + N^, (32) 

where the thermal contribution to the tadpole integral could be determined analytically at all temperatures because 
Aln = 0. The term Vq, as well as the vacuum contributions to the tadpole integrals (which depend on m a ) cancel 
when taking the difference. The critical temperature T c can be easily deduced from Eq. (|32p noting that 4>{T C ) = 0. 
The result is T c = y/12/N f n =V3f„. 



D. Nonlinear model in the chiral limit 



In the chiral limit of the nonlinear O(N) model both parameters e and h must be sent to zero. In the Hartree 
approximation and in the TR scheme, pions respect Goldstone's theorem by remaining massless in the phase of 
spontaneously broken chiral symmetry, see Fig. |6] In this phase, the cr-field is effectively frozen out due to its infinite 
mass. There is a first-order phase transition at a critical temperature T c = y3f n . At this temperature, the condensate 
drops to zero discontinuously, while the pion mass starts to increase continuously from zero above this temperature. 
In the restored phase, the a becomes degenerate in mass with the pions. This is the reason why T c assumes the same 
value as in the large- iV limit of the linear model. When inspecting the gap equations ([2"3")h we observe that they 
become identical with Eqs. (IM1) for M a = M„ = in the chiral limit and above T c (where <f> — 0). Therefore, we 
obtain the same equation (|32|) that determines the value of T c as in the linear case in the large- N limit. 

However, when applying the Hartree approximation in the CTR scheme no physical solutions can be obtained: the 
condensate goes to zero, <fi — > 0, and the masses of sigma and pion go to infinity, M CT , oo. This situation is 

similar to the nonlinear case with explicit symmetry breaking. 
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FIG. 6: The pion mass and the condensate as a function of T in the 0(4) nonlinear model in the chiral limit in the TR-scheme 
and Hartree-approximation for m CT — > oo (in practice m a — 250 MeV is used). The solid line corresponds to the physical case. 
The dashed line corresponds to the unstable or metastable solution of the gap equations. 
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FIG. 7: The pion mass and the condensate as a function of T in the 0(4) nonlinear model in the chiral limit using the 
LN-approximaion in the TR-scheme (full) and CTR-scheme (dashed) for m a — > oo. 



In the large- TV approximation, the phase transition is of second order with a critical temperature T c — VSfn, 
both in the TR and in the CTR scheme, see Fig. [7] Below T c the a mass is infinite, while the pions are massless, 
respecting Goldstone's theorem. Above the critical temperature the masses of the chiral partners become degenerate, 
M a = M n > in the TR scheme, and M a = M r = in the CTR scheme. At first sight, it is surprising that the a 
field becomes massless above T c . This behavior can be traced to our choice of the renormalization scale [i = m a — > oo. 
In fact, this is similar to what was observed in Ref. [l3[ (cf. Fig. 3 of that work), when increasing the renormalization 
scale in the large- TV limit in the CTR scheme. 

V. CONCLUSIONS 

In this work we have investigated the linear and the nonlinear O(N) model at nonzero temperature. An auxiliary 
field has been used to derive the effective potential. This method allowed us to establish a simple and mathematically 
rigorous relation between the linear and nonlinear versions of the model. This also leads to differences when comparing 
our results with previous treatments of the O(N) model, see below. The gap equations for the temperature-dependent 
masses and the condensate were derived using the C JT formalism. To regularize the divergent vacuum terms we applied 
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LN-CTR 


LN-CTR 


LN-TR 


LN-TR 


rriTT = m^ hvs 


m-x — s> + 




rrw ->■ 0+ 


m-x = m^ hvs 


nw — > + 


tutt = m£ ys 


m-x — > + 


lin 


second order at 
m CT ~ 500 MeV 




second order at 
m CT ~ 750 MeV 




crossover 


second order 
/To 


crossover 


second order 

/To" 
T -J—f 


nonlinear 


<S> 




first order 


first order 


no transition 


second order 


crossover 


second order 



TABLE I: The symbol (g) indicates that no reasonable result can be obtained due to tachyonic pion propagation. In those cases 
the phase transition becomes cross-over for smaller sigma masses and of first order for sigma masses higher than the shown 
values. m-x = m^ hva corresponds to the physical case of nonzero quark masses, m^ hvs — 139.5 MeV. 



the counter-term method (CTR scheme) as well as the so-called trivial regularization (TR scheme) where divergent 
terms are simply ignored. 

Table [J shows a compilation of the results for the various scenarios studied in this paper. The first row summarizes 
the results for the linear case, while the second those for the nonlinear case. In the first four columns we show the 
results for the Hartree approximation, the first two for the CTR and the next two for the TR scheme, for the case of 
explicit chiral symmetry breaking and in the chiral limit. The last four columns show the corresponding results for 
the large- TV limit. In the cases indicated with a <g>, we were not able to find physically acceptable solutions due to 
tachyonic pion propagation. In all other cases, we indicated the nature of the phase transition and, if independent of 
the <7 mass, the critical temperature. As one observes, T c = y/12/N = v^/w in the chiral limit for all scenarios, 
independent of the details (linear vs. nonlinear, or CTR vs. TR, or Hartree vs. large- N). In the cases where the order 
of the transition depends on the tr mass, we indicated the value of m a where the transition is of second order; it is 
crossover for smaller and of first order for larger values of m^. 

We now compare our results to previous work. In Ref. [Tj|, the O(N) model for N = 4 was studied in the CJT 
formalism without using the auxiliary field method. Although not studied in that work, we repeated the respective 
calculations varying the a mass. We find that, in the Hartree approximation and in the case of explicitly broken 
chiral symmetry, the phase transition changes from crossover to first order for m a ~ 940 MeV in the TR scheme and 
for ffifj ~ 680 MeV in the CTR scheme. This is consistent with our results obtained with the auxiliary field method, 
although the critical values for m a are somewhat larger for the method of Ref. [13j. In the chiral limit, the method 
of Ref. [l]| yields a first-order phase transition for all m a values. Furthermore, Goldstone's theorem is not fulfilled 
due to a non-vanishing pion mass in the phase of broken chiral symmetry. In the large- N limit, the results of Ref. 
fl3| coincide with ours, since the effective potentials are identical, cf. remark after Eq. (p28")l . 

The auxiliary field method has been applied previously to examine properties of the O(N) model to leading [27l.[28j 
and next-to-leading order in the 1/N expansion (2t| |30|. To leading order in 1/N the a and 7r fields have the same 
mass irrespective of whether chiral symmetry is explicitly or only spontaneously broken. Thus, in the chiral limit 
there are four instead of three massless bosons. The phase transition is of second order with a critical temperature 
of T c = y/U/Nf„. In the case of explicitly broken symmetry there is a crossover phase transition and four massive 
particles. To next-to-leading order including renormalization [29| the results change as follows: in the chiral limit 
there are three Goldstone bosons since the cr-field becomes massive. The phase transition is of second (or higher) 
order. In the weak-coupling limit the critical temperature is T c = y/12/ (N + 2)/ ff and above the critical temperature 
the masses of the chiral partners become degenerate. The key difference in our study to the afore mentioned Refs. 
[27L I28I ] is the correct treatment of the limiting process regarding the constraint imposed by the nonlinearity: the 
tr-mass is therefore infinite in the phase of broken symmetry. 

A natural next step is to include sunset-type diagrams in the CJT effective action, which lead to nonzero imaginary 
parts for the self-energy of the quasiparticles and, in turn, to a nonzero decay width [l6j |. Another project is the 
extension to nonzero chemical potentials [30j . A further interesting study would be the inclusion of additional scalar 
singlet states, since the nature of their constituency is quite unclear [34j |. Finally, the application of the auxiliary 
field method should also be instructive for more complicated systems incorporating additional vector and axial vector 
mesonic degrees of freedom 35] . 
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